Gravitational corrections to Yukawa systems 
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Abstract. We compute the gravitational corrections to the running of couplings in a scalar- 
fermion system, using the Wilsonian approach. Our discussion is relevant for symmetric as well as 
for broken scalar phases. We find that the Yukawa and quartic scalar couplings become irrelevant 
at the Gaussian fixed point. 
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I. INTRODUCTION 

The lack of renormalizability of Einstein's theory does 
not preclude the possibility of calculating quantum cor- 
rections to low energy processes due to graviton loops 
[l| . This effective field theory approach has been applied 
to calculate corrections to the gravitationalpotential |2:] 
and the running of Newton's constant [24a]- Graviton 
loops also contribute to the beta functions of matter cou- 
plings. This has been studied in the case of a scalar field 
in [6j . More recently, there has been considerable interest 
in (and controversy about) the corrections to the beta 
function of gauge couplings 0. Aside from the intrin- 
sic theoretical interest, such effects could have obvious 
applications to grand unified theories, whose characteris- 
tic energy scale is not too distant from the Planck scale, 
where gravity becomes strong. In fact, it has been argued 
recently [8| that in the determination of the GUT scale, 
quantum gravitational effects could be more important 
than two loop effects. 

With these motivations in mind, and in the same spirit, 
we will calculate here the gravitational effects on the beta 
functions of a simple Yukawa theory, consisting of one 
scalar and Nf fermion fields. We will do our calculations 
in flat Euclidean space, and therefore we will not calcu- 
late here the effect that the matter has on the running 
of the gravitational couplings (e.g. Newton's constant), 
but at least in the limit where the matter couplings are 
negligible, this effect is easily calculable |9|. 

In addition to the above, there is also another reason 
for studying this problem. If we look for a fundamental, 
as opposed to effective, theory of quantum gravity, there 
is now the concrete possibility that a purely field theo- 
retic solution can be obtained, provided that the renor- 
malization group has a fixed point with a finite number of 
UV attractive (relevant) directions. A theory with these 
properties is said to be asymptotically safe and has the 
same good properties (finiteness, predictivity) as, for ex- 
ample, QCD. The failure of perturbation theory means 
that the Gaussian fixed point of gravity does not have the 
desired properties. Work done in the last ten years has 
provided rather convincing evidence for the existence of a 
suitable nontrivial fixed point in pure gravity; see [lol| for 
reviews. It is then important to make sure that this fixed 
point persists also when interacting matter is brought in. 



In the case of scalar interactions, this was discussed in 
[Tl| . It was shown that there exists a "Gaussian mat- 
ter fixed point" , where the gravitational couplings are 
nonzero and slightly shifted relative to pure gravity, but 
all scalar selfinteractions are asymptotically free or zero. 
Our results imply that such a fixed point exists also in 
the presence of a Yukawa coupling. 

Finally we mention that asymptotic safety may play 
a role also in the standard model. Some evidence for a 
nontrivial fixed point in Yukawa systems has appeared 
recently If this was the case, then the calculations 
presented here are necessary to complete the picture by 
including also the gravitational interactions. 

II. YUKAWA SYSTEM 

In this section we set up the calculation. The flow of 
the renormalized couplings will be computed on a flat 
Euclidean background using an exact flow equation. An 
infrared cutoff, denoted k, is introduced via a cutoff term 
ASk, in order to define a scale dependent generating 
functional of connected Green's functions: 

-S[$]-/ J$~ASfc[$] 



M^fe[J] = -log /[d$]e 



(1) 



In flat space the cutoff term has the general form 
AS'fe[$] = y d^x<i>Rf{-d^)<i> and Rf{z) is constructed 
so as to suppress the contributions to the functional inte- 
gral from the infrared modes of the field $. For a scalar <j), 
we choose Rf.iz) — k'^r{z/k'^), with r{y) = {l — y)6(l — y) 
[l3| . leading to the substitution —d^ = z Pk{z) = 
z-\-k'^r{z/k'^), a kind of cutoff-propagator. For a fermion 
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The cutoff-corrected Legendre transform F^ = Wk — 
J d^xJcj)— ASk{4>) defines the effective average action Tk 
satisfying the renormalization group equation jl4i . Hsj 

' Rk] dtRk 



1 



dtTk = - STr 



(2) 



where t — Ink and STr denotes a functional trace, in- 
cluding a factor —1 for fcrmions. We will restrict our 
considerations to functionals F^ of the following form 

rfe[g^,.,0, ^, ^A] = Jd*x{Lb + Lf + Lg + LGF + Lgh). (3) 

The theory contains a single scalar field with Lagrangian 
Lb = Vg(iZ^V'^0V^0-f y((/.)) . 
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We choose the potential V to be even in 0. Then, 
there are Nf Dirac fermions ip with {/(iV^)-symnietric 
Lagrangian 

Lf = ^{^Z^i^rO^^^P - D^^pr^') + * H{q^) {p i,) . 
The covariant derivative is Df^ip = d^ij^ + ^u^cdJ'^'^'^', and 
D^'ip = dfj.'ip — ^ujficd'ipJ'^'^ , where w^cd is the spin connec- 
tion and J'^'^ — I [7'^, 7'*] are the 0(4) generators. We 
will choose the 0(4) gauge such that the vierbein is sym- 
metric, so that all vierbein fluctuations can be written 
in terms of the metric fluctuations and there are no 0(4) 
ghosts [1^. For the time being we keep the function H{(j)) 
general . On the other hand we will set = Z^ — 1 
and neglect anomalous dimensions. 

For gravity we have the Einstein-Hilbert Lagrangian 

where Z — l/(167rG). Similarly to previous analyses, 
we shall work with the background-field method. We ex- 



pand around constant backgrounds, which we still denote 
9 fill = S^}^, 0, Ip and V' with corresponding fluctuations 
/i^^, (p, X and X- For diffeomorphisms we fix a covariant 
background gauge, with gauge fixing term 
Z ^ „ 1+/3, 



L 



— AM'' pp. , ,, ,, ,, 



and the ghost action term consequently given by 



L 



gh 



We also employ the tensor decomposition 
^i^iy + 9f^v^ + di,v^ + {df^d^cr ~ jS^^d'^a) + \5^^h 
where d^'hj^^ = r|^'''hj^^ = 9^1^ = and /i = 5^'''h^,^, for 
tensor {hj^^), vector (w^) and scalar {a,h) fluctuations of 
the metric. 

The second order expansion of the Lagrangian ^ in 
the fluctuations /i^^^, w^, cr, h, ip, x and x is given by: 



2 \a I 16 



_3_ 
32 ~ 



-Zd"^ -2V - 2iH^^ aV + 3 



,/3 



-Zd''nd^h-—h 
16a 32 



/?2-3a 



Zd^ - 2V - 2iHijjijA h 



+iH'p{^Px + X^) + ^Hh{^x + Xyj) + ( ^d^vu 



3i 3z 

— a^9^cr - Y^d^h ] (iPl'^X - X7'''0) 
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where the primes denote derivatives w.r.t. (p and 7'^'^'' = 
{J^^jjP}. We also use the redefinitions 
> Wu, which remove the Ja- 



—d^a a and V—d^v^ — > 
cobians arising from the tensor decomposition. 

In order to write the RG flow for our sys- 
tem using ^ we introduce a supermultiplet T"^ ~ 
(/i^y , Vf^, Cfi,Cf^, cr, h, (fi, x'^ , x) containing all the field fiuc- 
tuations of the system and the matrix of operators 

Ffc^'' = sT^^'^Tr ■ "^^^ cutoff matrix Rk is chosen such 
that, adding it to FJ. leads to the replacement iS^ 



III. BETA FUNCTIONS. 

Let us define the dimensionless field (p — (p/k, and the 
dimensionless functions v{4>) = V{k(p)/k^ and h{(p) = 
H{k(p) /k. The running of V and H is obtained matching 
r - Jd'^x(^V + iHtP^P^ . Then, v = -iv + ^j)v' + k'^^V , 

and h = —h + cph' + k^^H . We present here the beta 
functionals for v and h, in the gauge P — 1 and expanding 
to first order in the dimensionless Newton constant G = 
k^G (the full expressions are nonpolynomial in G): 



V — —Av + t 



h = -h + (f)h' - 



Nf 



87r2 (1 + /i2) 



3 + 2v" 
327r2 {1 + v 

hh'^ (2 - 



_^i3-a)v'^2 + v") 



27r(l 



^v(3 + 2a)_ 



n{l+v"f 
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0{G'), 
hh'^ (4 



327r2 {l+v"y 

+ a(29 + 96/i2 + 48/i^) ~ , ,4a-6 - (3-2a)i;" + /i2(15-4a) + 2/i2(3-a) 



+ {2 + h^ 



((2 + /i2) «" + 2/i2A 



2n {1 + h^y {1 + v"Y 



+ 0{G'). (4) 



Fixing the form of the potentials and expanding around 
an appropriate basis of operators one may find the run- 



ning of any coupling of interest. We consider in the 
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following local power-law potentials, expanding either 
around ((/)) — or {(j)) = ^/k. Concerning h, from now 
on we restrict ourselves to a simple Yukawa interaction 
h = y0. 

a. Expansion around {(j)) — 0. For a quartic potential 



(5) 



inserting in Q we find, in the gauge a = and in the 
approximation Aq = 0, 



An = 



3 + 4A2 



Nt 



Ao = -2A 



A4 = 



327r2 (1 -I-2A2) 87r2 ' 

Nfy"" 3A4 



BGAo 



87r2 8^2(i + 2A2)2 ^(i + 2A2)' 



9A2 



Nfv' 



27r2(l + 2A2)^ 87r2 
+ 3GA 2 - 10^2 + 36A2 + 24A3 



0(G2), 



y 



^(l + 2A2)^ 



87r2(l + 2A2)' 



167r(l + 2A2 



In general, the beta functions would depend nonpolyno- 
mially on Aq and G. In the approximation Ao = 0, G 
appears only polynomially: the highest power of G oc- 
curs in A4 and is 2. In all other terms G appears at most 
linearly. 

When a 7^ one has to add the following correction 
terms: 

~ 29 + 180A2(l + A2) 
Ay = aGy- 



AA2 



AA4 



2aGA' 



2aGA4 



16^(l-f 2A2)' 
1-I-6A2(1 + A2) 
7r(l + 2A2)' 
1 + I4A2 



^(1 + 2A2)' 



(7) 



b. Expansion around a VEV. Depending on the sign 
of A2 , the potential ([5]) can be used to describe both the 
symmetric and the broken phase of the theory. In the 
latter case it may be more convenient to expand v around 
the VEV (0) = (k > 0), such that 



f'(V^) =0. 



(8) 



If we restrict ourselves to fourth order polynomials, v has 
the form 



+ 04(0' 



(9) 



The new couplings, in the broken phase where A2 < 0, 
are related to those in ([5]) by ^4 = A4, k = — A2/2A4, 6*0 = 
Ao — A2/4A4. The beta functions of these couplings can 
be derived from these relations and (|6]). Alternatively, 
one can obtain the running of n by deriving ([8]), which 
yields 



(10) 



For the broken phase, using Eq. ^ and retaining terms 
up to first order in G, we then obtain 

3 + I6K6I4 Nf 3G9o 



-A0O 
-2k ^ 



327r2 (1 + 8k6'4) 87r2 (1 + Ky^) 
3 Njy^ 



167r2 (1 

i2 



94k)2 16^2(1 + ^y2 



902 



Nfy^ 



27r2 (1 + 8k94 



y 



8772 (1 -h Ky2)^ 7r(l + 8K6l4)^ 

2 - I6K6I4 (3 + 8K6I4) 



167r2 (1 + ny^f (1 -|- 8k6'4)^ 
-3k?/2 (1 + 8K6I4 (7 + IQnOi)) - K^y^ (1 + 56K6I4) 



3Gy 



I66I4K (1 + 404k) 

16TT{l + y^KY{l + 894K)^^ 
-3y^K (1 + 86'4k) (9 + SO^k) + 192y^eiK^ (3 + 166'4k) 



+256y'^eiK^ (1 + 46'4k) 
We do not give here the 



(11) 

0{a) corrections to these for- 



mulae. We notice that unlike in the expansion around 
{(j)) — 0, here 9q appears only in its own beta function. 
Up to order G, there is no approximation involved in set- 
ting 00 = in the beta functions of k, 04 and y, as is 
natural in an expansion around flat space. 



IV. DISCUSSION. 

The standard MS result for the beta function of the 
Yukawa coupling is y = j^^pi + . . .. On the other hand, 
neglecting G and A2 in ^ or neglecting G and k in (fTTI) . 

we remain with y — + . . .. The difference is due to the 
fact that here we neglect the anomalous dimensions of (f) 
and ip. Since their contribution is not very small, our re- 
sults are not quantitatively accurate, but they should still 
give a reasonable qualitative picture of the gravitational 
corrections. We also stress that even though here we ana- 
lyze a toy model, our result for the leading one loop grav- 
itational correction applies also to realistic theories. In 
particular when the Yukawa couplings form a matrix y^ , 
every beta function i/ij will receive the same correction 
(27/167r)Gyij. The inclusion of anomalous dimensions 
is currently under study. Switching off the gravitational 
corrections, our results are in agreement with those of 
|12i | , when the anomalous dimensions are neglected. Fur- 
thermore, the results for Ai in ([S]) are also in agreement 
with those of ll']. We have given in Eqs. (|6|) and (fTTj) 
also the beta functions of the vacuum energy Ao and 0o. 
One can see the leading contributions, proportional to 
(3 — 4:Nf), the difference between the number of bosonic 
and fermionic degrees of freedom. 

Having used an expansion around flat space, gravity is 
off shell. This is the cause of the dependence of the results 
on the gauge parameter a (and /?, the dependence on 
which we have computed but not reported here for sim- 
plicity). We note that the sign of the leading corrections 
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does not change as long as a > 0; we have also checked 
that it remains the same at least for < /3 < 1.8, which 
comprises the most popular gauge choices. Furthermore, 
there are arguments showing that if a was allowed to 
run, a = would corresponds to a nonperturbative fixed 
point ,17]. This suggests that the results obtained for 
a = are probably the most reliable. 

The procedure also generically depends on the choice 
of cutoff scheme, and in particular on the cutoff func- 
tion r{y). The leading terms in the beta functions of 
A4 and y turn out to be independent on this choice, but 
not the gravitational corrections, which are related to a 
dimensionful coupling. In the results presented above 
we only used the cutoff r{y) = (1 — y)9{l ~ y), so the 
scheme dependence is not manifest, but the numerical 
coefficients of the gravitational correction would change 
if we used another cutoff function. We have checked that 
the leading gravitational correction is proportional to a 
single integral involving r{y), so that the ratio of the 
leading correction terms in ^ and (O is independent of 
r. Furthermore, the sign of the gravitational correction 
would be the same for any choice of r{y) that satisfies 
the boundary and monotonicity conditions to be a good 
cutoff. 

The system ([5]) has a (Gaussian) fixed point when A2 = 
X4 = y = 0. Without gravity both A4 and y are marginal, 
but the gravitational correction makes them irrelevant. 
In fact the critical exponents are 2 — (3 + 2a)G/Tr, 
— {3 + 2a)G/7T and — (27 + 29Q!)G/167r, corresponding to 
the eigenvectors A2 — 3A4/167r^, A4 and y respectively. 
(Note that the gravitational corrections depend on a but 
are always negative.) This is a remarkable result, be- 
cause in the standard model these couplings are free pa- 
rameters, to be determined by experiment, whereas here 
they are predicted to be zero at high energy. Any value 
they have at low energy is due to the nonlinearity of the 
RG flow. This result may change in the presence of other 
matter fields: it was shown in [llj that minimally coupled 
matter fields can change the sign of the critical exponent, 
making A4 relevant. Then its value at low energy would 
be a free parameter, while at high energy we would have 
asymptotic freedom. 

All this holds both for positive and negative A2. How- 
ever for negative A2 we may obtain an improved pertur- 
bation series [l^ expanding both v and h around the 
VEV. Then, the beta functions are those given in ([TT|). 



Most of the comments made above holds also in this 
case. The main difference lies in the fact that, in the 
absence of gravitational corrections, the fixed point now 
has 64 ~ y ~ and k — 3/327r^. Remarkably, the beta 
function of k does not receive any gravitational correc- 
tion, as was already noted in [6| for the potential (0) with 
^0 = 0, even taking into account the scalar field anoma- 
lous dimension. This is a general property: for any scalar 
potential v, using (fTO| and (|4|), 

o,, I V^^'" hNf^h' 

IGtt'^v" (1 + v"f 2(l + /i2)V2w" - ^' 

We stress again that the beta function of k obtained 
from the relation k — — A2/2A4 + A2A4/2A4 together with 
© has a G dependence in it. Also note that the general 
beta functional oih \xv ^ can be used to calculate the 
running of any term of the form 0"?/;?/;, in particular of 
an explicit fermionic mass. 

The gravitational corrections are of order G — 
fc2/Mpjjjjjj.jj and therefore can be treated perturbatively 
at low energies. They may not be negligible at the GUT 
scale, though. Beyond the Planck scale the gravitational 
corrections seem to be large and unbounded. The theory 
may still be meaningful provided all couplings (in par- 
ticular G) reach a fixed point. It is known that in the 
Einstein-Hilbert truncation gravity has a nontrivial fixed 
point, also in the presence of minimally coupled matter 
fields. Since the Yukawa system has a Gaussian fixed 
point, one can conclude that the theory of gravity cou- 
pled to scalars and fermions also has a fixed point, which 
we may call a "Gaussian matter" fixed point. However, 
it is clear that to study the properties of this fixed point, 
in particular the critical exponents, it is necessary to cal- 
culate also the beta function of G. There is also the pos- 
sibility that the matter sector exhibits a nontrivial fixed 
point [12]. Preliminary results indicate that, as long as 
G, < 1, this fixed point would also exist in the presence 
of gravity. We plan to discuss these matters in more 
detail elsewhere. 
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